A theory of magnetization reversal in nanowires 



Robert S. Maier 

Depts. of Mathematics and Physics, University of Arizona, Tucson, AZ 85721, USA 

ABSTRACT 

Magnetization reversal in a ferromagnetic nanowire which is much narrower than the exchange length is believed 
to be accomplished through the thermally activated growth of a spatially localized nucleus, which initially oc- 
cupies a small fraction of the total volume. To date, the most detailed theoretical treatments of reversal as a 
field-induced but noise-activated process have focused on the case of a very long ferromagnetic nanowire, i.e., 
a highly elongated cylindrical particle, and have yielded a reversal rate per unit length, due to an underlying 
assumption that the nucleus may form anywhere along the wire. But in a bounded-length (though long) cylin- 
drical particle with flat ends, it is energetically favored for nucleation to begin at either end. We indicate how 
to compute analytically the energy of the critical nucleus associated with either end, i.e., the activation barrier 
to magnetization reversal, which governs the reversal rate in the low-temperature (Kramers) limit. Our treat- 
ment employs elliptic functions, and is partly analytic rather than numerical. We also comment on the Kramers 
prefactor, which for this reversal pathway does not scale linearly as the particle length increases, and tends to a 
constant in the low-temperature limit. 

Keywords: Kramers theory, weak spatiotemporal noise, Landau-Lifshitz-Gilbert equation, activation barrier, 
Kramers prefactor, ferromagnetic cylinder, ferromagnetic nanowire, magnetization reversal, Neel-Brown theory 

1. INTRODUCTION 

The phenomenon of magnetization reversal is of intense interest to the data storage and recording industries, 
for obvious reasons. Most treatments of this phenomenon have employed the standard techniques of micromag- 
netic modeling: analytically or numerically solving a partial differential equation satisfied by the magnetization 
vector, which is valid in an averaged, coarse-grained sense, on scales large compared to the crystal lattice pQ. 
Ferromagnetic samples of many geometries have been studied theoretically [2]. 

Provided that a sufficiently strong reversed field is applied to a ferromagnetic sample, the reversal of the 
direction of magnetization is an essentially deterministic process, which if the sample is sufficiently large (larger 
than the coherence length) will begin with the formation of a spatially localized 'nucleus' of reversed magneti- 
zation. This nucleus will grow, eventually yielding a global reversal of the magnetization. This is in contrast to 
the case when the applied field is less strong than the nucleation field, in which a critical nucleus can only be 
formed by thermal fluctuations, though its subsequent growth may be deterministic. The two types of reversal 
are sometimes called field-activated and noise-activated. In many applications, such as quantifying the likelihood 
of loss of stored data, it is the second phenomenon, more closely tied to statistical mechanics than to the theory 
of deterministic evolution equations, which is of interest. 

The difference between the two types of reversal is made clear by an analogy. Consider an oversimplified 
overdamped scalar model of the form rii = —U'{m) + noise, where U(-) is a bistable effective potential, the two 
minima of which correspond to an alignment of the magnetization with and against an applied field. As the 
field strengthens, this double- well potential may become sufficiently skewed that the latter of the two minima 
disappears, leading to a rapid fall into the remaining minimum, which corresponds to alignment with the field. 
This is the essence of field-activated nucleation. On the other hand, if the field is insufficiently strong, then 
there will remain two minima: a metastable one (against the field) and a stable one (with the field). It is 
clear that noise-induced nucleation is in essence a standard problem of noise-induced escape from a metastable 
state. Of course, a satisfactory model will need to be more sophisticated than this scalar model. Even if the 
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sample is small ("zero-dimensional"), so that reversal takes place by a process of spatially coherent rotation of 
the magnetization vector, the vectorial nature of the magnetization is crucial. 

The process of noise-activated reversal by coherent rotation was first treated theoretically by Brown [3], who 
by adapting the Kramers theory of weak-noise escape rates was able to work out an asymptotic (low-temperature) 
expression for the reversal rate in a single-domain particle with uniaxial anisotropy. His work has been extended to 
the case of cubic anisotropy. However, the issue of which reversal mechanisms are dominant in large ferromagnetic 
particles remains controversial. The simplest sort of non-zero-dimensional sample is a ferromagnetic nanowire: 
a long, thin cylinder of width at most a few tens of nanometers. The most extensive treatment of field-induced but 
noise-activated magnetization reversal in such a sample is that of Braun 4, 5 6 , |S] . By applying the infinite- 
dimensional Kramers theory of Langer and Coleman j^j to a stochastically perturbed Landau-Lifshitz-Gilbert 
equation governing the evolution of the magnetization vector along the nanowire, he worked out analytically an 
asymptotic expression for the weak-noise reversal rate. However, his treatment assumes a toroidal sample (with 
periodic boundary conditions), and applies only in the limit of infinite sample length. In this limit, the nucleus 
of his model is a metastable field configuration of soliton type, with a simple functional form. It comprises a 
'critical droplet' of reversed magnetization bounded by a tt wall and a —tt wall; this droplet may form with 
equal likelihood anywhere along the sample. His computed reversal rate is accordingly proportional to the 
sample length L. It also contains a factor exp(— AE/ksT), where AE is the energy of formation of the droplet. 
Like Brown's asymptotic reversal rate formula, Braun's contains corrections to the Arrhenius (i.e., exponential) 
low-temperature falloff. For reasons explained below, both fall off as T -1 / 2 exp(— AE/ksT) as T — > 0. A non- 
Arrhenius behavior of the reversal rate as T — > has been found in stochastic simulations of single-domain 
particles, at least [Til] . 

Aharoni [111 I12| has criticized Braun's not including the (nonlocal) magnetostatic self-energy in his Landau- 
Lifshitz-Gilbert Hamiltonian. If the magnetization is not solenoidal (V ■ m ^ 0), there will be a magnetostatic 
volume charge, which will generate a demagnetizing field. This term can be neglected in a planar topology |13|. 
but not in the cylindrical topology. Aharoni has also criticized the assumption that the magnetization vector 
is constant across the cross-section of the cylindrical sample, which prevents reversal from occurring via a 
'curling' mode. However, this assumption is presumably justified in the limit of small sample width. Also, 
by what amounts to dimensional analysis, Braun has shown that in this limit the nonlocal part of the self- 
energy is relatively small, and the local part may be absorbed into a renormalization of crystalline anisotropy 
constants |14| . With these emendations, Braun's analysis seems to be applicable to thin cylindrical samples of 
at least some materials, though of course the assumption of a very long sample with periodic boundary conditions 
is unrealistic and will need to be dropped. 

In this extended abstract we indicate how the analytic treatment of Braun can be extended to the case of 
a thin cylinder with flat ends, at which open (i.e., Neumann) boundary conditions on the magnetization vector 
are appropriate. We do this by adapting our work on reversal in a real Ginzburg-Landau equation perturbed 
by weak spatiotemporal noise |15| . Since we do not take the L — > oo limit, our expression for the metastable 
nucleus configuration is more complicated than Braun's. It involves elliptic functions [lfi| . which have long been 
used in the analytic modeling of periodic micromagnetic structures. (See ref. 1171 §2.4, and ref.^3]) Sophisticated 
numerical simulations at the magnetic moment level of noise- activated reversal, in a thin cylinder with flat ends, 
show that a droplet of reversed magnetization forms at one of the ends and spreads inward |19II20| . This reversal 
mode is in agreement with our analytic treatment. We also comment on the prefactor in the low-temperature 
reversal rate, i.e., the factor that multiplies exp(— AE /ksT). As L — + oo, it does not grow linearly in L, but 
rather tends to a constant, since the formation of a droplet of reversed magnetization away from both ends of 
the cylinder is energetically disfavored. We also find that the prefactor includes no T~ x / 2 factor. That is, the 
rate of noise-activated magnetization reversal, if physically reasonable boundary conditions are imposed, follows 
an Arrhenius law. 

2. ESCAPE FROM A METASTABLE STATE 

The theory of noise-activated fluctuations away from a metastable state [21], including escape to a stable state, 
was developed most significantly by Kramers |22| . Major contributions were also made by Eyring, who worked 
on escape in multidimensional state spaces in the context of reaction-rate theory. It should be emphasized that 
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the Kramers approach to escape rates is asymptotic: increasingly accurate in the weak-noise limit, in which 
escape per unit time becomes exponentially unlikely. In this limit, the phenomenology of noise-activated escape 
from the domain of attraction of a metastable state involves eventual motion along an 'optimal trajectory' that 
typically goes uphill from the metastable state to a saddle point, near which exit occurs. By definition, this is 
the escape pathway of least resistance. The escape rate is computed as the integral of the outgoing probability 
flux through a region on the boundary of the domain of attraction which is centered on the saddle point. The 
optimal trajectory picture and the associated computations become increasingly valid in the weak-noise limit. 

The present author and his collaborators have extended Kramers theory to the case when the system being 
modeled, if noise is absent, has 'nongradient' dynamics: for example, when its order parameter x (not a mag- 
netization, in general) is finite-dimensional but does not evolve according to any law of the form x = — VE/(x). 
This could be because the dynamics include dissipation or non-conservative forces. In the nongradient case 
the phenomenology of escape may be quite different from that described above. The domain of attraction of 
a metastable state may have an unstable limit cycle as its boundary, and the outgoing optimal trajectory may 
spiral into this limit cycle. This is in fact a feature of the abovementioned coherent rotation model of Brown [3], 
which is Hamiltonian but includes a damping term. In two-dimensional models in which the boundary of the 
domain of attraction is an unstable limit cycle, the Kramers escape rate must be computed as the integral of the 
outgoing probability flux over the entire boundary, rather than merely over the vicinity of a saddle point (which 
may not be present). This usually gives rise to anomalous weak-noise behavior of the Kramers prefactor, i.e., 
the factor that multiplies the usual exponential suppression factor in the asymptotic weak-noise escape rate. 

In this section we review the extensions to Kramers rate theory that we shall need in understanding and 
extending the results of Braun 0J |SJ El [7| |S] , with an emphasis on the prefactor and on the phenomenology of 
escape. We shall use dimensionless units throughout. Suppose a point x in finite-dimensional Euclidean space 
(or possibly on a surface such as a sphere) evolves according to the overdamped Langevin equation 

x = u(x) + e 1 / 2 ^)- (1) 
Here u is a drift field, rj is normalized white noise, which satisfies 

(»fc(ti)»&(t2)} = M(*i-*a), (2) 

and e <C 1 is a noise-strength parameter. Suppose the drift field has a metastable point S and a stable point S', 
with a saddle point denoted U (for unstable) between them, on the common boundary of their domains of 
attraction. If x = S initially, the time for noise-activated escape can be defined as the time needed for x to 
reach some specified region R that lies strictly within the domain of attraction of S', and contains S' itself. This 
escape time t ox n(S — > S') is a random variable, the distribution of which in the weak-noise limit will not be much 
affected by the choice of region R. As e — > 0, its mean (t cxit (S — > S')) will grow exponentially (to leading order), 
and its distribution will increasingly resemble that of an exponential random variable. The escape rate Ts~*s' 
can be defined as the reciprocal (t CX i t (S — > S")) -1 . Just as there is noise-activated motion from S to <S", there 
is noise-activated motion from S' back to S, but reverse motion is exponentially suppressed in the weak-noise 
limit, relatively as well as absolutely. 

The escape rate Tg^s' may be written as r° exp(— AW/e), where AW is the minimum amount of energy 
the noise must put into the system to drive it over the saddle U. By definition, the second factor is the Kramers 
factor, and F° is the Kramers prefactor. The case when the drift field u is the negative gradient of a potential 
function, i.e., u = — VV^, is the easiest case to treat, since the optimal escape trajectory is simply a time-reversed 
relaxational trajectory, and satisfies x = — u(x). It turns out that W = 2V (the factor of 2 is sometimes absorbed 
into a redefinition of e) . The weak- noise behavior of the escape rate can readily be expressed in terms of linear 
approximations to the drift field at x = S and x = U, by the Kramers-Eyring formula. In the two-dimensional 
case, this formula is 



IW ~ ^ \/\ A|| (U) | A,, (S) J ±±L± exp(-2A V/ e ) , e^O. (3) 

Here A||(5) = —dus.\\/dxs,\\ is the eigenvalue of the linearized negative drift field at S whose correspond- 
ing eigenvector points along xgii, the direction locally parallel to the departing optimal trajectory. Similarly, 
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X±(U) = —duu,\./dxu,±. is the eigenvalue of the linearized negative drift field at U whose corresponding eigen- 
vector points along Xj/Xj the direction locally perpendicular to the approaching optimal trajectory, and so forth. 
A||(Z7) is the eigenvalue of the linearized negative drift field at U that corresponds to the unstable, or 'downhill' 
direction; it is negative. The formula generalizes immediately to the case when the order parameter x has 
arbitrary finite dimensionality d. In general, the prefactor will contain a product of d — 1 factors resembling 

\f \±(U) ' one ^ or eacn direction transverse to the optimal trajectory. 

The standard formula contains a constant prefactor, which does not include any power of the noise 
strength e. But the derivation of © makes clear how such a power could arise. The stationary probability 
density function (PDF) of the order parameter x is po(x) = exp(— W(x)/e), up to normalization. The slowest 
decaying nonstationary mode of the drift-diffusion equation for the PDF is a quasistationary PDF pi, which 
characterizes how probability equilibrates between the domains of attraction of S and S'. Like po, it will contain 
an exp (— 2V(x)/e) factor, so it will typically be concentrated on the (^(e 1 / 2 ) transverse lengthscale near the saddle 
point U. But unlike po, p\ will have opposite signs in the two domains. Along the x^n direction through U, it 
will change sign. The rate Ts~,s' is the rate at which probability diffuses across the saddle, from S toward S'. As 
Kramers showed, this is easily evaluated as an integral of the probability flux. In this integral there are several 
powers of e that must be kept track of. For example, there is a normalization factor proportional to e~ d / 2 in pi 
(as there is in po)- Also, it follows from the Langevin equation JQ) that besides drifting along u, the PDF diffuses 
with diffusivity e. A careful evaluation of the probability flux integral yields the Kramers prefactor shown in J2J, 
with all powers of e canceling. 

A prefactor containing a nontrivial power of e could arise from the potential energy V = V(x) behaving non- 
quadratically near x = S or U, causing one or more of the four eigenvalues appearing in @ to equal zero |23j. 
As mentioned, non-Arrhenius behavior will also arise in the case of a two-dimensional order parameter if the 
separatrix between the domains of attraction of S and 5" contains no saddle point, and is an unstable limit 
cycle. For this to be the case, the drift u cannot be the negative gradient of a potential function. A transverse 
integration of the probability flux generated by the quasistationary PDF p\ over an O(l) lengthscale around 
a limit cycle, rather than over an 0(e 1/ ' 2 ) transverse lengthscale near a saddle point, will yield an escape rate 
larger by a factor roughly proportional to e 1 / 2 . In such a situation, 

r s _s' ~ const x e" 1 / 2 exp(-AW/e) e -> , (4) 

where 'const' can be evaluated as an integral around the separatrix, will replace J3J). 

All of the preceding remarks apply to the case of finite-dimensional overdamped dynamics, where the Langevin 
equation for the order parameter has the form l[T)). If there is damping but the strength of the damping is not 
taken to infinity, a more appropriate Langevin equation in dimensionless units is 

x + ax = u(x) + e 1/2 77(f), (5) 

where a > is a damping constant, or equivalently 

v = -av + u(x) +e 1 / 2 r 1 (t), ^ 

where v = x. This is effectively an overdamped Langevin equation for a 2d-dimensional order parameter (x, v) 
that diffuses only in the velocity direction. But if a has some fixed positive value, the weak-noise behavior of the 
Kramers prefactor in this non-overdamped model is surprisingly similar to that discussed above. The prefactor 
becomes a-dependent [3T], but there is still no e-dependence if u = -W with V quadratic near S and U. 
When d = 2, an e -1 / 2 factor will appear in the Kramers prefactor if u is not derived from a potential and there 
is no saddle on the separatrix, for essentially the reason already discussed. There are clearly some subtleties 
associated with the reduction of the non-overdamped model of (JSJ to the overdamped model of QJ when a — > oo, 
but we do not discuss them here. 

Infinite-dimensional Kramers theory is a theory of rate and relaxation phenomena exhibited by random fields, 
rather than by randomly evolving finite-dimensional order parameters. However, it is formally quite similar to 
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finite-dimensional Kramers theory. This is perhaps not entirely clear from standard expositions UJ, but we shall 
make it so. For simplicity, consider a field <f> — 4>{x, t) perturbed by noise of strength e, which evolves with 
overdamped dynamics on a spatial domain [0, L) according to 

:j) = -SH/S(/> + e^ 2 r]{x,t). (7) 

The first term on the right-hand side is the functional derivative of an energy functional, and rj(x,t) is normalized 
spatiotemporal white noise, which satisfies 

(rj(xi,ti)ri(x2,t2)) = 5(xi - x 2 )5{t 1 - t 2 ), (8) 

i.e., has a flat power spectrum in space as well as time. Typically 

< L V(*) 2 



n\d>] = 
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v{4>{x)) 



dx, (9) 
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so that the evolution equation is a stochastically perturbed Ginzburg-Landau equation, namely 

4> = Ac/) - V'{(j>) + e 1/2 r)(x, t). (10) 

If the potential V is bistable or multistable then noise-induced transitions between stable spatially extended 
states, each of which is concentrated around a field value that is a minimum of V, may occasionally take place. 
Each such transition will involve a passage through an intermediate (unstable) saddle state. 

These spatially extended states are strongly affected by the boundary conditions at x = and x = L. 
Suppose that V{<j>) has two local minima, at cj> = S, S 1 , with the latter being the global minimizer; and that 
the intervening local maximum is at cj> — U. If periodic or Neumann boundary conditions are imposed, one 
would expect that the metastable, saddle, and stable states are the constant field configurations cfi = S, U, S"; 
though if (say) Dirichlet boundary conditions are imposed, one would expect them to differ, especially near x = 
and x = L. Actually, even in the periodic and Neumann cases the saddle state most relevant to noise- induced 
transitions may not be the constant solution cf> = U. If L is sufficiently large, there may be a lower-energy saddle 
in the field configuration space, which is spatially oscillatory rather than constant ^H]- The importance of such 
nonconstant saddle configurations was first pointed out by Langer |24| in a study of the L — > oo limit of escape 
from a metastable state. Due to their oscillatory character they are sometimes called "Langer's bounces" . Each 
such oscillatory saddle configuration can be viewed as a background metastable (4> = S) state, within which a 
droplet of the stable (cf> = 5") state has formed. The droplet is critical, which means that if pushed in the right 
direction, the saddle configuration may roll 'downhill' to the stable rather than the original metastable state. 

The Kramers-Eyring formula generalizes readily to the infinite-dimensional case, becoming 



r~^J "r'fii ex P (-2AH/e), e^O. (11) 

Here {A„((/> S )}^L and {A„((/>„)}^L are the eigenvalues of the linearizations of the negative drift field 8TL/5<j) 
at the metastable state 4> = <f> s (x) and the dominant saddle <f> = 4>u(x), arranged in increasing order. The only 
negative eigenvalue in either set is Xq((/) u ), which corresponds to 'downhill' motion away from the saddle. It is 
easy to see that in the two-dimensional case, ifTTJl reduces to ©• 

Because the numerator and denominator of the quotient inside the square root are products of an infinite 
number of eigenvalues, which tend to infinity as n — > oo, each diverges. However, their quotient can be defined 
as a limit of finite-dimensional truncations. Usually it is recast as a determinant quotient, computed from the 
linearized deterministic evolution operators at (f> = <p s , <f> u |251 126j . Consider for example a small perturbation £ 
away from the metastable state, i.e., <fi — <fi s +£. Then to leading order £ = — A[<^ s ]£, where A[cj) s ] = S 2 Tt/6(t) 2 [(j) s ] 
specifies the linearized zero- noise dynamics near <fi = <p s . Similarly, A[<^„] specifies the dynamics near cf> = <p u . 
For the model defined by ©, A [t/> a ], A [</>„] equal -A + V" {4> s {x)) , -A + V" (<t> u {x))- One can write 




r*.^ . ~ ' T \ , ; t , exp(-2AH/e), e . (12) 
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The two determinants of Schrodinger operators are not denned individually, but their quotient can be made 
sense of by any of several regularization schemes; for example, zeta- function and dimensional regularization |26|. 

If the boundary conditions at x = 0, L are periodic, so that the spatial domain [0, L] is effectively a circle, 
and the dominant (lowest-energy) saddle state <fi = (j> u {x) is not the constant solution cf> = U but rather a field 
configuration containing a critical droplet, the asymptotic formula Ijl2(l must be modified. Since the droplet 
configuration can be shifted arbitrarily around [0,L], the operator A[<p u ] will have a zero eigenvalue \i(<fi u ), the 
accompanying eigenfunction being a Goldstone mode that expresses translation invariance. So in this case the 
prefactor of i|12|) will diverge, signalling a change in the phenomenology of escape. In the infinite-dimensional 
space of field configurations, the separatrix between the domains of attraction of the metastable and stable 
states will now include an unstable limit cycle of length L, and noise-activated escape may take place anywhere 
along it. We have already explained how to handle such separatrices. Since the escape rate must be computed 
by integrating the outgoing probability flux around the separatrix, rather than on a transverse lengthscale of 
magnitude 0(e 1 ^ 2 ), the prefactor will be larger by a factor that grows like Le" 1 / 2 as L — » oo. The asymptotic 
weak-noise behavior of the escape rate in this case, when L is fixed, is accordingly of the form 

r^ a , ~ C(i) e - 1 /2 e xp(-2(AH)(i)/6), e^O, 
C(L) ~ const x L, L — > oo , 

where 2(A7i)(L) is the energy barrier, i.e., the energy of formation of the critical droplet within the spatial 
domain [0, L\. This will tend to a limiting value as L — > oo, since the droplet tends to a limiting shape. We stress 
that the asymptotic linear growth of the prefactor with L in this case occurs for a simple physical reason: 
if periodic boundary conditions are imposed, the critical droplet may form with equal likelihood anywhere in the 
spatial domain. Departures from linearity will occur when L is small, but they are due to the droplet assuming 
its limiting shape only in the large- L limit. The large- L behavior of the prefactor displayed in is well-known, 
but the finite- volume case of fixed L has only recently been systematically treated ^3 • 



3. STOCHASTIC MICROMAGNETISM 

We can now relate the low-temperature micromagnetic reversal rate theories of Brown [3] and Braun ^|S1|SJ[7||S] 
to finite-dimensional and infinite-dimensional Kramers theory, respectively, and explain the physical origin of the 
T -1 / 2 prefactor that modifies each of their Arrhenius laws. Our discussion will indicate how Braun's treatment 
of a very long (L — > oo) ferromagnetic nanowire can be modified and extended to the case of a thin ferromagnetic 
cylinder with flat ends, and fixed length L. 

The starting point for continuum micromagnetic modeling 1 is the Landau Lifshitz -Gilbert equation for 
the bulk magnetization vector m. If dimensionless units are used and the magnetization unit is chosen so that 
the saturation magnetization of the ferromagnetic material being modeled equals unity, i.e., so that |m| = 1 at 
each point in the sample, this equation is 

rh = —7m x h c ff + am x m, (14) 

where 7 > is the gyromagnetic ratio. The first term causes m to precess about the local magnetic field h c ff. 
The second is a phenomenological damping term introduced by Gilbert |2Zj , which causes m to relax toward the 
direction of h c ff . This becomes clear if the equation is rewritten as 

(1 + a 2 )7 -1 m = m x h c ff — am x (m x h e ff), (15) 

which follows by evaluating the cross product of m with each term. The Gilbert damping term, which is 
now universally used as a replacement for a less satisfactory damping term proposed earlier by Landau and 
Lifshitz [2HI j is still slightly controversial, in part because of the difficulty of calculating the material-dependent 
dimensionless damping coefficient a > from first principles. 

The formula for h g- must capture much of the physics. If the dynamics are essentially Hamiltonian (with 
the exception of the Gilbert damping), then h c ff = — SH[m]/Sm for some Hamiltonian functional H [m] of the 
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magnetic field configuration. The Hamiltonian employed by Braun for a thin cylindrical particle whose axis 
extends from x = to x = L, with m taken to depend only on x, is 

H[m] = f { [(d x m x ) 2 + {d x m v f + [d x m z f] + (m, Km) - h cxt • m} da. (16) 
Jo 

The first term is the exchange term, of quantum-mechanical origin, which favors the alignment of nearby magnetic 
moments, and the second is a quadratic form on the sphere S 2 3 m that models the anisotropy of the material, 
if any. For convenience the unit of length here is chosen to equal the exchange length, so the coefficient of the 
exchange term is unity. 

This Hamiltonian does not include a term of a type arising from classical electromagnetism, the omission 
of which has been criticized by Aharoni ^2 ^] . If V • m ^ in the interior of the sample then there will be 
a magnetostatic volume charge, which will generate its own magnetic field, with which the bulk magnetization 
will interact. The additional field is called the 'demagnetizing' field, for an obvious reason. In general, this 
m-generated field will add a term to the energy density which is nonlocal in space. However, as mentioned in 
Section ^ Braun |14j has given convincing arguments for replacing this term in the thin-cylinder limit by a local 
term, which can be absorbed into a renormalization of the anisotropy tensor K. 

Braun considered the effects on the Landau-Lifshitz-Gilbert equation H15|l of adding a small random term of 
the form e 1//2 rj(x,<) to h c ff, where J](x,t) is vectorial spatiotemporal white noise. Before discussing his results, 
it is useful to review Brown's much earlier theoretical treatment of ferromagnetic reversal in a single-domain 
particle. The equation H15|l makes sense as the evolution equation of a single vector m, regarded as the spatially 
uniform saturation magnetization across such a particle. The appropriate Hamiltonian is not (|16l) but rather 

H[m] = (m, Km) - h cx t • m , (17) 

so 

h cff = -5H[m}/5m = h oxt - 2Km . (18) 

The anisotropy tensor K should include both shape and crystalline anisotropies. Brown considered the case of 
uniaxial anisotropy, with the external field directed along the 'easy axis' of the particle. If this easy axis points 
in the x-direction, K = diag(— K e , 0, 0) for some K e > 0, so let h cxt = (h cxt ,0,0) with h cxt > 0. Then as a 
function of the angle 9 between m and the x-axis, 

H[m] = -K e cos 2 9 - h cxt cos 9. (19) 

Provided that h cxt < 2K e , this function on the sphere S 2 3 m has two minima: the metastable state 9 = it 
(antiparallel to h cxt ) and the stable state 9 — (parallel to h cxt ), which has lower energy. Their energies are 
respectively — K e ± h cxt . 

It is easy to check that according to the evolution equation 115f) . the particle magnetization m will spiral 
into one or the other of these two states. On the sphere, the separatrix between their domains of attraction 
is the parallel of latitude singled out by the equation 9 = cos _1 (— h cxt /2K e ). This is an unstable limit cycle. 
If a temporal noise term of the form e 1 / 2 r7(t), with e <C 1, is added to h c ff , then the noise will eventually drive 
the magnetization vector away from the metastable state and across this separatrix. It can be shown that the 
optimal trajectory for noise-activated escape is, in fact, a time-reversed infalling spiral. So the mechanism of 
activated magnetization reversal in this model is of the sort discussed in Section [21 However, it should be noted 
that when /i ext > 2K e , the metastable state becomes unstable, and there is no separatrix to be crossed. Reversal 
in the presence of a sufficiently high external field is field-activated rather than noise-activated. 

By examination, the 'height' of the unstable limit cycle in terms of energy is h 2 ^ t /AK e , so the energy barrier 
to reversal is 2/S.TL = 2 [h 2 xt /4K e — (—K e + /i cx t)] ■ The weak-noise reversal rate will be given by (@J, i.e., 

r - const x e- 1/2 exp(-2A?i:/e), e^0. (20) 

The factor 'const' can be calculated by integrating along the limit cycle the rate at which the probability density 
of m is absorbed. However, this derivation already makes it clear why the Nccl-Brown reversal rate, as it is now 
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called, falls off in the limit of low temperature in a non-Arrhenius way. The anomalous power e" 1 / 2 , i.e., T" 1 / 2 , 
arises from the need to integrate a probability flux around the entire limit cycle, rather than over an O^e 1 ^ 2 ) 
transverse lengthscale near a saddle point, as is more typical in the Kramers theory of noise-activated escape. 
This derivation can be extended to the case of a single-domain ferromagnetic particle with mildly non-uniaxial 
anisotropy tensor K, in which, as here, the barrier to reversal is an unstable limit cycle. 

We now turn to the thin nanowire computations of Braun. In their full generality, they apply to a thin 
nanowire of a ferromagnetic material with cubic anisotropy, the easy axis of which is aligned with the axis of the 
wire; which may be taken to be the cc-axis. By adding an appropriate multiple of the identity matrix to K (which 
has no physical effect, since |m| = 1 at all points in the sample), and rotating the coordinate system about the 
wire, K can be taken to equal diag (— K e , 0, Kh) for some K e ,Kh > 0. The z-axis will be a hard axis, and the 
x—y plane an easy plane, within which the magnetization vector at every point along the wire is preferentially 
confined. A hard axis of this sort could arise from a shape anisotropy rather than a crystalline anisotropy, e.g., 
from the cross section of the wire being elliptical rather than circular. 

For simplicity we consider here only the limiting case Kh — > oo, in which all fluctuational motion of the 
magnetization vector out of the x—y plane is suppressed. (Taking out-of-easy-plane fluctuations into account 
would alter the Kramers prefactor by multiplying it by an extra determinant quotient, but would not affect the 
qualitative nature of its e — > and L — > oo asymptotics.) For this case the effective Hamiltonian is 



H[m] 



{ [(d x m x ) 2 + (d x m y ) 2 ] - K e m 2 - h cxt m x } dx. 



As a functional of the angle 9 between m (in the x—y plane) and the x-axis, 



nm 



— I — K e cos 2 9 — hext cos 9 



dx. 



(21) 



(22) 



In the absence of thermal noise, the stationary configurations of the deterministic Landau-Lifshitz-Giibert 
equation (|15fl will be those for which 8H [m]/<5m = 0, i.e., SH[9]/S9 = 0. That is, they will be those for which 
the magnetization inclination angle 9 as a function of distance x along the wire satisfies 



d 2 9 

2—r + 2K e sin 9 cos t 

dx z 



h ext sin 9 = 0. 



(23) 



For the moment, suppose that the wire can be viewed in some sense as toroidal, i.e., suppose that periodic 
boundary conditions 9(0) = 9(L), 9'(0) — 9'(L) can be imposed. Then the solutions of iffify will include 9 = it 
and 9 = 0, which are the familiar states of metastable and stable magnetization, with the magnetization pointing 
against or along the external field, which is itself applied parallel to the wire. The solution 9 = cos~ 1 (—h cxt /2K e ) 
is also stationary in time, though unstable. This is the uniform (Neel-Brown) barrier to magnetization reversal, 
which in this context is a saddle point in the infinite-dimensional space of field configurations. 

However, if L is sufficiently large, there will also be a nonuniform barrier to magnetization reversal. By 
proceeding as in our analysis |15| of the Ginzburg-Landau equation (|10l) , one can show that 



9 u {x) 



a\ sn 2 ((x — xo)/c\m) + ciq 
bi sn 2 ( (x - x )/c \ m) + b 



(24) 



where do, ffli, bo, b%, c, m depend on L and h ext /2K e , and xq is arbitrary, is also an unstable stationary solution, 
and is at lower energy than the Neel-Brown barrier. Here sn(-|m) is the Jacobi elliptic function with modular 
parameter m € (0, 1), the period of which is 4K(m), where K(m) is the first elliptic integral ^B]. The parame- 
ter m is of course chosen so that 9 U , as a function of x, has period L. If such a lower-energy barrier to reversal 
exists, its 'height' in terms of energy, which may be written as ATt(L), will determine the exponential falloff of 
the rate of escape from the metastable state 9 = tt, as the noise strength tends to zero. 

At this point, we have gone significantly beyond the analysis of Braun. It is not difficult to check that the 
solution 9 = 9 u {x) can be viewed as a critical droplet: for many choices of parameter, it resembles a droplet of the 
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9 = phase inserted into a background consisting of the 9 — n phase, and delimited by a ir wall and a — n wall. 
However, its shape depends on the length L of the wire. As L — > oo, it tends to a limiting shape, which may 
be expressed in terms of hyperbolic trigonometric functions. This is possible because as L — > oo, the modular 
parameter m tends to unity, and sn(x|m) degenerates to tanhir. The formula (|24() is new in this context, though 
the limiting shape as L — > oo was obtained by Braun. Similar expressions for periodic micromagnetic structures, 
involving elliptic functions, have been derived in other contexts |29| . 

It follows immediately from the discussion in Section that for this model of a thin cylindrical particle of 
length L with periodic boundary conditions, the magnetization reversal rate T satisfies 

r ~ C(L) c-i/2 exp(-2(AW)(L)/e), e - , 
C(L) ~ const x L, L — » oo . 

The factor e 1 / 2 is due to the fact that the separatrix between the metastable and stable domains of attraction 
contains a limit cycle of length L, since the droplet may be formed anywhere along the wire, on account of 
translation invariance. (xq above is arbitrary.) This is also responsible for the asymptotic linear growth of the 
prefactor with L. These two features were noted by Braun, but we now see they have a simple physical origin. 

The present approach, in which the shape of the critical droplet for fixed L is precisely specified in terms 
of elliptic functions, yields precise values for the exponent 2AH(L) and the prefactor C(L). Details of the 
calculations, which are numerical rather than analytic, may appear elsewhere. What is of greater interest, 
however, is the modification of this analysis to include more realistic boundary conditions. Braun has pointed out 
that in a cylinder of finite length, 'open' boundary conditions with the magnetization perpendicular to the ends 
may be appropriate [H) . In the present scheme this corresponds to taking 0'(O) = 0'(L) = 0, i.e., to Neumann 
boundary conditions. 

Equipped with Neumann boundary conditions, the Euler-Lagrange equation (|23[1 still has the metastable 
and stable states 9 = n, 9 = as solutions. A detailed analysis, which will be published elsewhere, reveals that 
by choosing parameters appropriately, the unstable elliptic solution 6 = u (x) of Q24JI can also be made to satisfy 
Neumann boundary conditions. However, when this is done, it becomes a 'half-droplet' delimited by a ir wall 
and either x = or x = L. Unlike the case of periodic boundary conditions, nucleation at low temperatures is 
therefore not likely to take place anywhere along the wire, but must, on energetic grounds, preferentially begin 
at either end. Once a critical half-droplet has formed near either end, it may roll 'downhill' in the magnetization 
configuration space by spreading along the wire, replacing the original value 9 — ir by 9 — as it goes. 

The important point is that if Neumann rather than periodic boundary conditions are imposed, there will 
be no infinite degeneracy of the transition state, no unstable limit cycle, no non-Arrhenius 0(e -1 / 2 ) behavior of 
the Kramers prefactor, and also no asymptotic linear dependence of the prefactor on L as L — > oo. The noise- 
activated magnetization reversal of a long thin ferromagnetic cylinder actually fits much better into traditional 
Kramers theory if Neumann boundary conditions are imposed. 

ACKNOWLEDGMENTS 

This research was supported in part by National Science Foundation Grant No. PHY-0099484. 

REFERENCES 

[1] J. Fidler and T. Schrefl, "Micromagnetic modeling — the current state of the art," J. Phys. D 33, pp. R135- 
R156, 2000. 

[2] A. Aharoni, Introduction to the Theory of Ferromagnetism, Oxford University Press, Oxford, 2nd ed., 2000. 

[3] W. F. Brown, Jr., "Thermal fluctuations of a single-domain particle," Phys. Rev. 130(5), pp. 1677-1686, 
1963. 

[4] H.-B. Braun, "Thermally activated magnetization reversal in elongated ferromagnetic particles," Phys. Rev. 
Lett. 71(21), pp. 3557-3560, 1993. 



9 



[5] H.-B. Braun, "Fluctuations and instabilities of ferromagnetic domain-wall pairs in an external magnetic 
field," Phys. Rev. B 50(22), pp. 16485-16500, 1994. 

H.-B. Braun, "Statistical mechanics of nonuniform magnetization reversal," Phys. Rev. B 50(22), pp. 16501- 
16521, 1994. 

H.-B. Braun and H. N. Bertram, "Nonuniform switching of single domain particles at finite temperatures," 
J. Appl. Phys. 75(9), pp. 4609-4616, 1994. 

[8] H.-B. Braun, "Kramers' rate theory, broken symmetries, and magnetization reversal," J. Appl. Phys. 76(10), 
pp. 6310-6315, 1994. 

[9] L. S. Schulman, Techniques and Applications of Path Integration, ch. 29. Wiley, New York, 1981. 

[10] E. D. Boerner and H. N. Bertram, "Non-Arrhenius behavior in single domain particles," IEEE Trans. 
Magnetics 34(4), pp. 1678-1680, 1998. 

[11] A. Aharoni, "Magnetostatic energy of the soliton in a cylinder," J. Magnetism and Magnetic Materials 
140-144, pp. 1819-1820, 1995. 

[12] A. Aharoni, "Comment on "Kramers' rate theory, broken symmetries, and magnetization reversal" ," J. Appl. 
Phys. 80(5), pp. 3133-3134, 1996. 

[13] J. S. Broz, H.-B. Braun, O. Brodbeck, W. Baltensperger, and J. S. Helman, "Nucleation of magnetization 
reversal via creation of pairs of Bloch walls," Phys. Rev. Lett. 65(6), pp. 787-789, 1990. 

[14] H.-B. Braun, "Nucleation in ferromagnetic nanowires — magnetostatics and topology," J. Appl. Phys. 85(8), 
pp. 6172-6174, 1999. 

[15] R. S. Maier and D. L. Stein, "Droplet nucleation and domain wall motion in a bounded interval," Phys. 



Rev. Lett. 87(27), paper no. 270601, 2001. Available as |arXiv:cond-mat /0108217 



[16] M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, Dover, New York, 1965. 

[17] W. F. Brown, Jr., Micromagnetics , no. 18 in Interscience Tracts on Physics and Astronomy, Interscience, 
New York, 1963. 

[18] A. Aharoni, "Reduction in coercive force caused by a certain type of imperfection," Phys. Rev. 119(1), 
pp. 127-131, 1960. 

[19] G. Brown, M. A. Novotny, and P. A. Rikvold, "Micromagnetic simulations of thermally activated mag- 
netization reversal of nanoscale magnets," J. Appl. Phys. 87(9), pp. 4792-4794, 2000. Available as 
arXiv:cond-mat/9909136 



[20] G. Brown, M. A. Novotny, and P. A. Rikvold, "Langevin simulation of thermally activated magnetization re- 
versal in nanoscale pillars," Phys. Rev. B 64, paper no. 134422, 2001. Available as arXiv:cond-mat/0101477 

[21] P. Hanggi, P. Talkner, and M. Borkovec, "Reaction-rate theory: Fifty years after Kramers," Rev. Modern 
Phys. 62(2), pp. 251-341, 1990. 

[22] H. A. Kramers, "Brownian motion on a field of force and the diffusion model of chemical reactions," Physica 
7(4), pp. 284-304, 1940. 

[23] C. Kipnis and C. M. Newman, "The metastable behavior of infrequently observed, weakly random, one- 
dimensional diffusion processes," SIAM J. Appl. Math. 45(6), pp. 972-982, 1985. 

[24] J. S. Langer, "Theory of the condensation point," Ann. Physics 41(1), pp. 108-157, 1967. 

[25] J. S. Langer, "Statistical theory of the decay of metastable states," Ann. Physics 54, pp. 258-275, 1969. 



10 



[26] A. J. McKane and M. B. Tarlic, "Rcgularization of functional determinants using boundary conditions, 
J. Phys. A 28(23), pp. 6931-6942, 1995. Available as |arXiv:cond-mat / 9509126 



[27] T. L. Gilbert and J. M. Kelly, "Anomalous rotational damping in ferromagnetic sheets," in Proceedings 
of the 1st Magnetism and Magnetic Materials Conference, AIEE Special Publications (T-78), pp. 253-263, 
(Pittsburgh), 1955. 

[28] S. Iida, "The difference between Gilbert's and Landau-Lifshitz's equations," J. Phys. Chem. Solids 24, 
pp. 625-630, 1963. 

[29] M. J. Shirobokov, "On the mechanism of cobalt magnetization," Dokl. Akad. Nauk SSSR 24(5), pp. 426-429, 
1939. 



11 



